Computation of dynamical correlation functions of Heisenberg chains in a field 
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We compute the momentum- and frequency-dependent longitudinal spin structure factor for the 
one-dimensional spin-1/2 XX Z Heisenberg spin chain in a magnetic field, using exact determinant 
representations for form factors on the lattice. Multiparticle contributions are computed numerically 
throughout the Brillouin zone, yielding saturation of the sum rule to high precision. 



Exact solutions of quantum models, either on the lat- 
tice or in the continuum, are invaluable in the study 
of nonperturbative aspects of low-dimensional physics. 
Bethe's construction of the full set of eigenstates of 
Heisenberg's spin exchange modeli, by the method today 
commonly known as the Bethe Ansata^, has ultimately 
led to a wide variety of exact results on the thermody- 
namics of integrable models, with a broad spectrum of 
applications^. 

A long-standing limitation of the Bethe Ansatz was 
that the dynamics of such models was not directly ac- 
cessible. In recent years, however, enormous progress 
has been made for correlation functions of the Heisen- 
berg spin chain 5 ' 6 ' 7 . In particular, matrix elements of 
any local operator between two Bethe states can now be 
written as matrix determinants 6 -. Combined with formu- 
las for eigenstate normal, this yields exact expressions 
for form factors on the lattice, thereby opening the door 
to the computation of dynamical correlation functions. 

In this paper, we implement this program for one of 
the cornerstones of the theory of integrable models, the 
anisotropic Heisenberg antiferromagnetic chain in a mag- 
netic field: 
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with periodic boundary conditions. The anisotropy pa- 
rameter A will be chosen to lie in the gapless regime 
— 1 < A < 1, and therefore the model provides a well- 
controlled realization of quantum critical behaviour. 

Our interest lies in the space- and time-dependent spin 
correlation functions. We will numerically compute the 
longitudinal dynamical spin structure factor, which is de- 
fined as the Fourier transform of the spin-spin correlation 
function: 
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where the subscript c means that we take the connected 
part. This quantity is directly accessible experimentally 
through neutron scattering (see e.g. 0^3)- 

The longitudinal and transverse correlations for the 
isotropic Heisenberg chain in a field were studied at fixed 
momentum q in and the transverse ones for XX Z 



at zero field and at q = n in [T^]- More recently, the 
two-particle contributions to the longitudinal structure 
factor for the XX Z chain in a field at q = ir/2 were stud- 
ied for all energies in [l^. In particular, this allowed (at 
low energies) for a numerical check of conformal scaling 
exponents^ to an impressive degree of precision. Here, 
we present results for all momenta q, including multipar- 
ticle contributions. This yields data beyond the reach 
of conformal field theory, and allows us to quantitatively 
evaluate the precision of our results via the sum rule. 

The exact solution through the Bethe Ansatz for model 
JTJ is well-known 2 . The reference state is taken to be the 
state with all spins up, |0) = T)i- The Hilbcrt 

space separates into subspaces of fixed magnetization, 
determined from the number of reversed spins M. We 
take the number of sites N and of reversed spins M to 
both be even. Eigenstates in each subspace are com- 
pletely determined for 2M < N by a set of rapidities 
{Xj}, j — 1, M, solution to the Bethe equations 
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where A = cos £. Each choice of a set of distinct half- 
integers {Ij}, j = 1,...,M (with Ij defined mod(TV)) 
uniquely specifies a set of rapidities, and therefore an 
eigenstate. The energy and momentum of a state are 
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The ground state is given by I 6 

In terms of form factors, the structure factor can be 
written as a sum 



S"(g, W ) = 27rX;i<0|S| 

a=i0 



q)\ 2 5{uj - uj a ) 



(5) 



over the whole set of eigenstates (distinct from the 
ground state) in the fixed M subspace. Each term in 
© can be obtained 6 as a product of determinants of M- 
dimensional matrices, fully determined for a given eigen- 
state by a knowledge of the corresponding set of rapidi- 
ties. For the sake of brevity we do not reproduce the 
expressions for these matrices here. 



2 





0.5 1 1.5 2 2.5 0.5 1 1.5 2 2.5 







A = 0.3 




2p 











0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 
CO 





FIG. 1: Constant momentum q slices of the structure factor as a function of frequency ui, at q = tt/2 and Sn/4, for XXZ 
with A = 0.3 and A = 0.7 at M = N/4. 2p denotes two-particle contributions, while 4p denotes four-particle contributions. 
Resulting curves are extrapolations to infinite size from N = 4096 (2p) and 384 (2p + 4p). 2p contributions carry most of the 
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FIG. 2: Longitudinal structure factor as a function of momentum q and frequency u, for XXZ with A = 0.3, and A = 0.7 
at M — 7V/4. Two-, four- and six-particle contributions (labeled 2p, 4p and 6p) are plotted separately, for system sizes 
TV = 768, 128 and 80 respectively. 



Obtaining the spin structure factor thus involves three 
steps: scanning through the eigenstates, solving the 
Bethe equations, and computing the determinants to ob- 
tain the form factor. The fact that actually permits one 
to achieve a high degree of precision, is that the contri- 
butions are rapidly decreasing functions of the number 
of particles involved. Contributions from unbound states 



to the structure factor come from even numbers of one- 
particle excitations. Two-particle contributions are those 
for which there is only one hole in the distribution of I's 
in the interval corresponding to the ground state. Simi- 
larly, four- and six-particle contributions are those with 
two and three holes in the distribution of the I's. It is 
well-known that the Bethe equations © yield, for certain 
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TABLE I: Contributions of two-, four- and six-particle sectors 
to the sum rule for A = 0.3 (%). 



N 


2p 


4p 


6p 


Total 


768 


70.29 
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99.11 
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82.17 


17.52 


0.23 


99.92 



TABLE II: Contributions of two-, four- and six-particle sec- 
tors to the sum rule for A = 0.7 (%). 
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Total 


768 


69.85 
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choices of Fs, complex rapidities in the form of strings 
(bound states). This leads to reductions in the deter- 
minants (we will publish the formulas elsewhere) . In the 
presence of a high field, the bound-state contributions are 
strongly suppressed. We compute the two-string contri- 
butions to show that this is indeed the case. The zero- 
field case, where bound states carry more weight, will be 
studied in a separate publication. 

We present results for two different cases: a chain with 
A = 0.3, and another with A = 0.7, both in a field such 
that M = N/A. For each case, we give two slices of the 
structure factor at fixed momentum q = 7r/2, 37r/4 as a 
function of frequency in Fig. ^ including contributions 
from up to four particles (six-particle and bound-state 
contributions are too small). We subsequently give den- 
sity plots for the whole region q 6 [0, n] in Fig. [21 for the 
same cases, including up to six-particle contributions. In 
all cases, we use eq. (JSJ where we broaden the delta 
functions to width log N/N in order to obtain smooth 
curves. The isotropic case A = 1 again requires special 
treatment, which we will present elsewhere. 



Computing the structure factor for all momenta allows 
us to check the saturation of the sum rule 

LXJ ^ L J 

For M = N/A, the sum rule requires C = 0.1875. 

We compute two-particle contributions to the struc- 
ture factor for system sizes up to N = 768. Four- and 
six-particle ones are computed up to N = 128 and 80 
respectively. Two-particle contributions are monotoni- 
cally decreasing with system size, whereas four- and six- 
particle ones first increase, then decrease, due to the in- 
terplay between the finite-size gap and state counting. 
The contribution to the sum rule from states with one 
bound state and two holes is computed up to N = 256, 
and yields only 4.4 x 10~ 7 % for A = 0.3, and 3.0 x 10~ 7 % 
for A = 0.7. Relative contributions are presented in Ta- 
bles H] [n] Two-particle terms account for about 70% of 
the sum rule, close to the 72.89% result for the isotropic 
chain in zero fielcU^i^. Including four-particle contribu- 
tions saturates around 99% of the sum rule for 128 sites. 
Our results for 80 sites, including up to six-particle con- 
tributions, saturate well over 99% of the total sum rule. 

In summary, we have numerically obtained correlation 
functions of integrable models on the lattice using deter- 
minant representations for form factors, in the fundamen- 
tal case of the anisotropic Heisenberg chain in a field. We 
have shown that by including multiparticle contributions, 
the sum rule can be saturated to very high precision. 
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